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1 Introduction 



Suppose that during a de Sitter stage of expansion the coupling constant of an Abelian 
gauge field evolves in time. Thus the kinetic term of the gauge field can be written (in 
four space-time dimensions) as 

S cm = -\j d 4 xV^Gf(<P)F a(3 F^, (1.1) 

where G is the determinant of the space-time metric, is a scalar field (which can depend 
upon space and time) and g(<p) = /(0)~ 1,/2 is the coupling |. This type of vertex is typical 
of scalar-tensor theories of gravity |IJ and of the low energy string effective action ||. 
Early suggestions that the Abelian gauge coupling may change over cosmological times 
were originally made by Dirac [[| (and subsequently discussed in J3J and in 0) mainly in 
the framework of ordinary electromagnetism. 

The dynamics of the field will be described by the action of a minimally coupled 
(massive) scalar. If the field is displaced from the minimum of its potential during in- 
flation, there will be a phase where the field relaxes. Provided the scalar mass is much 
smaller than the curvature scale during inflation such a phase could be rather long. During 
the de Sitter stage, the specific form of the expanding background will dictate, through 
the equations of motion, the rate of suppression of the amplitude of <fi. 

While the field relaxes toward the minimum of its potential, energy is pumped from the 
homogeneous mode of to the gauge field fluctuations. The function /(0) can be either 
an increasing function of (leading to a decreasing coupling) or a decreasing function of 
(leading to an increasing coupling). In both cases, depending upon the parameters of the 
model, the two-point correlation function of magnetic inhomogeneities increases during 
the inflationary stage. This implies that large scale magnetic fields can be potentially 
generated. 

The implications of large scale magnetic fields for cosmology, astrophysics and high- 
energy physics have been recently reviewed in || where various useful references are 
reported. The interested student will find in a morte physical introduction to the 
subject of large scale magnetic fields. The toy model discussed in the present paper is 
only an useful exercise suggesting a possible connection between the variation of gauge 
couplings and the production of large scale magnetic fields. 

Recently || the variation in the fine structure constant has been suggested by observa- 
tions of absorption lines at different frequencies. The results of these claims can be summa- 
rized by saying that a ern was smaller in the past. Defining, indeed Aa em = a to d ay — a P ast> 
measurements suggest that 

= (0.72±0.18)1(T 5 . (1.2) 

Q^em 

2 The Heaviside electromagnetic system of units will be used throughout the investigation. The effective 
"electron" charge will then be given, in the present context, by e(</>) = ei/(0) -1 / 2 . 
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To avoid confusions, in the present scenario, gauge couplings are constant today and their 
variation occurs prior to big-bang nucleosynthesis (BBN). 

The gauge coupling has to be (almost) constant by the time when the Universe is 
approximately old of one second, namely by the time of BBN. The abundances of light 
elements are very sensitive to any departure from the standard cosmological model. Hence, 
fluctuations in the baryon to photon ratio, matter-antimatter domains, anisotropies in 
the expansion of the four space-time dimensions, can all be successfully constrained by 
demanding that the abundances of the light elements are correctly reproduced. Following 
the same logic the variation in the gauge couplings can also be constrained from BBN 

The plan of the present contribution is then the following. In Section II the basic ideas 
concerning the model of evolution of the gauge coupling will be introduced. In Section 
III bounds coming both from the homogeneous and from the inhomogeneous evolution of 
4> will be described. In Section IV the evolution of the magnetic inhomogeneities will be 
addressed along the various stages of the model with particular attention to the role of the 
two-point function. In Section V the large scale magnetic fields produced in the scenario 
will be estimated. Section VI contains some concluding remarks. The considerations 
reported in the present paper extend and complement the results of [11]. 

Finally, in concluding this Introduction, I wish to remind a couple of useful references 
of two other lecturers of the Chalonge school. In Ref. an interesting discussion on the 
possible ambiguities arising in the observations of large scale fields in clusters is reported. 
In H another idea on the generation of magnetic fields is reported. 



2 Basic Equations 

Thanks to the high degree of isotropy and homogeneity of the observed Universe, the 
background geometry can be described using a (conformally flat) Friedmann-Robertson- 
Walker (FRW) line element 

ds 2 = G^dx^dx" = a 2 (r])[dr] 2 - dx 2 }, (2.1) 

where rj is the conformal time coordinate and G^ v is the four-dimensional space-time 
metric. The cosmic time coordinate (often employed in this investigation) is related to r] 
as a(rj) di] = dt. 

From the anisotropies of the Cosmic Microwave Background (CMB) it is consistent 
to assume that the Universe underwent a period of inflationary expansion of de Sitter or 
quasi-de Sitter type. Therefore a(r/) ~ —Vi/v during a phase stopping, approximately, 
when the curvature scale was Hi < 10 _6 Mp. For t] > r]i (possibly after a transient period) 
the Universe gets dominated by radiation [i.e. a(r]) ~ 77] and then, after decoupling, by 
dust matter [i.e. a(r]) ~ rf]. 
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The action describing the dynamics of the (Abelian) gauge coupling in a given back- 
ground geometry can be parametrized as 



S = Jd 4 x V^G 



,<j)-V{(j>)- 



\ J? 

V(0) . 



(2.2) 



In Eq. (|2.2|) <p is not the inflaton but an extra (scalar) degree of freedom taking care of 
the evolution of the gauge coupling. 

From Eq. ( |2.2j ) the equations of motion can be derived 



1 



0, 



V^GG^d u 
F 



9 2 {<t>) 



+ 



0. 



dV 



1 dg 



d<p 2^ 3 (0) d<p 



(2.3) 
(2.4) 



If the gauge coupling does not change, the evolution of Abelian gauge fields is con- 
formally invariant. Hence, using the the conformal time coordinate, the appropriately 
rescaled gauge field amplitudes obey a set of equations which is exactly the one they 
would obey Minkowski space. In order to simplify the explicit form of the equations of 
motion the electric and magnetic fields will be rescaled in such a way that the obtained 
system of equations will reproduce the usual (conformally invariant) system in the limit 
g — > constant. The rescalings in the fields are 



B = a 2 B, E = a 2 S, A = aA, 



(2.5) 



where B, £, A, are the flat-space quantities whereas B, E, A, are the curved-space ones. 
The explicit form of our system becomes, in the metric ( |2.1| ): 



d^± + m d± + dV a2 
drf di] dip 

dB 
dr] 




dg 



2g 3 a 2 



E 2 -B' 



— = -VxE, 



-E 



dr] lg 2 ((j)) 

V • B = 0, 

2dg_ 
gd<f> 



V -E 



9 2 W 



E ■ V0, 



V x B 



2dg 



V<p x B 



(2.6) 

(2.7) 

(2.8) 

(2.9) 
(2.10) 



where v is the bulk velocity of the plasma. The spatial gradients used in Eqs. ( |2.6| )-( |2~T(]| ) 
are defined according to the metric ( [2.1| ). In Eq. ( |2.11| ) the quantity 7i = din a/ drj has 
also been introduced. 7i is the Hubble factor in conformal time which is related to the 
Hubble factor in cosmic time as 7i = H/a where H = dlna/dt. 
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Once the background geometry is specified we are interested in the situation when the 
gauge field background is vanishing and the only fluctuations are the ones associate with 
the vacuum state of the Abelian gauge fields. Hence, Eqs. (|2.6Q - (|2.10|) allow to compute 
the evolution of <fi and the associated evolution of the two-point function of the gauge 
field fluctuations. 

Suppose that is originally displaced from the minimum of its potential. As far as 
the zero mode of is concerned the system of equations can be further simplified: 



d 2 



dV 



+ m-£- + —of 

or] 

= -Vxi 



drf " dr] ' d(f) 
dB 
drj 
d_ 
dr] 



l_dg_ 



2g 3 a 



E 2 — B 2 



9 2 (<P) 



E 



9 2 (<P) 



V x B, 



V • B = 0, V • E = 0, 



(2.11) 

(2.12) 

(2.13) 
(2.14) 



By now combining together the modified Maxwell's equations we obtain the evolution of 
the magnetic fields 

B" - 2—B' - V 2 B = 0, (2.15) 
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where the prime denoted derivation with respect to the conformal time coordinate (the 
over-dot will denote, instead, derivation with respect to cosmic time). 

Once the evolution of the metric is specified, Eq. ( |2.11| ) dictates a specific evolution for 
(f) and the evolution of (f) will determine, in its turn, the evolution of the gauge fields. The 
interesting initial conditions for the system are the ones where the classical gauge field 
background vanishes. Thus, when the homogeneous component of starts its evolution 
during the de Sitter phase, quantum mechanical fluctuations will be postulated as initial 
conditions of gauge inhomogeneities. 

When the background geometry evolves from the de Sitter phase to the subsequent 
epoch, massive quanta of are produced. The amount of the produced inhomogeneous 
modes of can be computed and it will be shown that the associated energy density will 
always be smaller than the one of the homogeneous mode. This analysis will be one of 
the subjects discussed in Section III. 



3 Constraints on the evolution of the gauge coupling 

Sub-millimiter tests of the Newton's law show that no deviations are observed down 
to distances as small as 0.1 mm fll^f . Therefore, should be massive. Of course the 
potential of <fi may be much more complicated than the one provided by a simple mass 
term. However, for sake of simplicity, a massive scalar will be analyzed since, already 
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in this case, interesting effects can be analyzed. In spite of the fact that this choice is 
apparently simple, various constraints on the scalar mass appear. 



3.1 Evolution of the homogeneous mode 

Suppose that the potential term driving the evolution of the gauge coupling is simply 



v{4>) 



m 



During the inflationary stage of expansion the scale factor evolves as a(r]) = 
i] < — 771. The evolution of is obtained by solving 

0" + V + ^0 = 0, 

where /1 = mjH\ and where the relation 7i ~ rj^ 1 ~ aH has been used. If // -C 1, for 
f] < —r]i the solution of Eq. ( p.2|) can be written as 



(3.1) 
-r)i/r)) for 

(3.2) 



Hv) =<f>i-fa[ — 



(3.3) 



The end of the inflationary stage of expansion may not be directly followed by the radi- 
ation dominated phase. In the intermediate phase the scalar mass is still small than the 
curvature scale but the curvature decreases, in general, faster than during the inflation- 
ary phase since the background is neither of de Sitter nor of quasi-de Sitter type. The 
evolution of the scale factor can be parametrized as a{rf) ~ r] a where, in order to fix the 
ideas, a ~ 2 could be assumed f\ The evolution of will simply be 



4>ih(v) = 01 + 02 



2a -4 



+ 



2a — 1 2a — 1 vr/i 



V \ 



2a- 1 



Vl < V < Vr 



(3.4) 



where the continuity between Eq. (|3.3|) and Eq. (|3.4| ) has been required, so that <j>i(—r)x) ■ 

(j) A (vi) an d <&(-Vi) = 0rhO/i)- 

After i] r the background enters a radiation dominated phase and the evolution of 

can be explicitly solved in cosmic time. The equation for , in this phase, is given by 



+ 3#0 + m 2 

3 

which in terms $ = a 2 0, becomes 

$ + nr 



0, H = -, 

a 



3 -H - V 

2 4 



$ = 0. 



(3.5) 



(3.6) 



3 The case a — 2 corresponds to a matter-dominated intermediate stage. 
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In the radiation-dominated stage of expansion Eq. ( ^6|) becomes 



m 2 + 



16t 2 



$ = 0, 



(3.7) 



whose solution can be written in terms of Bessel functions 



$(mt) 



AYi(mt) + BJi(mt) 



(3.f 



For mt <C 1, has a constant mode and a solution as t -1 / 2 . Recalling the relation 
between cosmic and conformal time and imposing the continuity of <fi and 
the solution of Eq. (|3.4|) the following form can be obtained: 



(m r} r 



with 



>1 + 



2a -4 6(1 



2a 



2a 



2a-l 



2a-l 



(3.9) 



which is valid for r] r < r] < i] m . The time rj m marks the moment where H ~ m. When 
mt > 1, the regime of coherent oscillations takes over and the solution (|3.8|) implies that 



the energy density stored in <fi decreases as a -3 , meaning that (j) c (v) ~ rf z l 2 . Since the 
coherent oscillations decrease as a~ 3 there will be a typical curvature scale H c and a 
typical time r] c at which the coherent oscillations become dominant with respect to the 
radiation background. This moment is determined by demanding that 



H 2 M 2 ( - 



2/2 f a m \ ' 



(3.10) 



which also implies that 

H c ~rri(p 4 , (3.11) 

where (p = (pi/Mp. Eq. (|3.11|) has been obtained without tuning the asymptotic value of 
<fi to the minimum of its potential. If such a tuning is made, the amplitude of oscillations 
at r] m will be smaller than 0i and it will be given, according to Eq. (|3.9| ), by 02(^i/^r) 2a_1 - 
Thus, the scale H c will be defined by a different relation namely: 



2 j.2 

m 02 



El 

#1 



2(2n-l) 
a + 1 



m 

W r 



leading, ultimately, to 



rn 



H 2 M 2 



4(2a-l) 



m 



(3.12) 



(3.13) 



In the approximation of instantaneous reheating [i.e. r] r ~ 771] , H r ~ ffi. Therefore, 



from Eq. (|3.13|) H c is smaller than the value determined in Eq. (|3.11|) by a factor 



(m/Hi) 2 . In the approximation of matter-dominated reheating (i.e. a ~ 2), the result 
of instantaneous reheating is further suppressed by a factor (H T /HiY as one can easily 
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argue from Eq. (|3.13|) . From Eqs. ( |3.5|) -( [3~6"1) , the evolution of will go as r\ 3 when 



coherent oscillations start dominating. 

In spite of the possible tunings made in the asymptotic values of 0, after r) c there will 
be a typical time at which the field will decay. In order not to spoil the light elements 
produced at the epoch of BBN has to decay at a scale larger than H ns ~ T 2 s /Mp (where 
T ns ~ MeV. Since is only coupled gravitationally the typical decay scale will be given 
by comparing the rate with the curvature scale giving that 



3 

Tfl 

H ^ T ~M2> h ^ ( 3 - 14 ) 

implying that m > 10 4 GeV. This requirement also demands that the reheating temper- 
ature associated with the decay of will be larger than the BBN temperature. 

In order to illustrate some concrete examples of the various possibilities implied by our 
considerations, suppose that m ~ 10 3 TeV and suppose that the asymptotic value of 
is fine-tuned to its minimum. Furthermore, suppose that the reheating is instantaneous. 
Then, according to the picture which has been presented, inflation stops at a scale Hi ~ 
10 13 GeV and starts oscillating at a curvature scale H m ~ 10 3 TeV. The coherent 
oscillations will then become dominant at a curvature scale H c ~ 10~ 8 GeV (having 
assumed 02 ~ Mj>). The coherent oscillations of will last down to H^ ~ 10~ 20 GeV. 
After this moment the Universe will be dominated by the radiation produced in the decay 
of 0. Notice, for comparison, that the BBN curvature scale is H ns ~ 10 -25 GeV so that 
the decay occurs well before BBN (five orders of magnitude in curvature scale). 

Another illustrative example is the one where m ~ 10 6 TeV. In this case the decay of 
occurs prior to the EWPT epoch, namely 

> H ew . (3.15) 

In fact H ew = v^T^/Mp ~ KT 17 GeV (with N cS = 106.75 and T ew ~ 100 GeV) 
whereas, from Eq. ( |3.14| ), ~ 10~ 9 GeV. In more general terms we can say that in 
order to have the decay occurring prior to the EWPT epoch we have to demand that 
> H ew which means that m > 10 5 TeV. 

In closing this section two general comments are in order. If no fine-tuning is made in 
the asymptotic amplitude of 0, the typical scale of the coherent oscillations will almost 
coincide with m. However, the possibility ip 1 is still left if, for some reason, we want 
H c < m. 

The decay of and the consequent freezing of the gauge coupling should occur prior 
to the EWPT epoch and the baryon number should be generated, in the present context, 
at the electroweak time. Suppose, for example, that this is not the case and that the BAU 
has been created prior to the electroweak scale. Suppose, moreover, that the decay of 
occurs after baryogenesis. Then the temperature of the radiation gas before the decay of 
will be ~ T m (a m /a0) ~ m(m/M P ) 5 . Thus, the entropy increase due to the decay of 
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(f) will be AS ~ (Tdocay/T^) 3 where Td ocay ~ H^Mp. This implies that AS ~ m/Mp. It 
has been observed in different contexts that in order to preserve a pre-existing BAU one 
should have AS < 10 5 [ffSJ, [T(|. Thus, this bound would imply m > 10 14 GeV. This is 
the reason why the present analysis will assume that the decay of occurs prior to the 
electroweak time and that the BAU is generated at the EWPT or shortly after. 



3.2 Evolution of the inhomogeneous modes 

When the Universe passes from the inflationary stage to the subsequent radiation dom- 
inated expansion, inhomogeneities of the field <fi are generated. This may invalidate the 
original assumptions and introduce further complications by adding qualitatively new 
constraints on the scenario. 

It is useful to recall that the inhomogeneities of (f) can be interpreted, in the framework 
of second quantization, as quanta of the field 0. Hence, the inhomogeneities produced 
because of the sudden change of the geometry from the de Sitter epoch to the radiation 
dominated epoch can be counted by estimating the number of quanta produced by the 
sudden change of the geometry according to the well known techniques of curved space- 
times Wh - 



Consider the first order fluctuations of the field <fi 

(p(x,r}) = <p(rj) + 5(j)(x,r]), (3.16) 

whose evolution equation is, in Fourier space, 

if>" + 2Hip'+ [A; 2 + mV]^ = 0, (3.17) 

where ip(k, rj) is the Fourier component of 5<j)(x, rf). In order to count the number of quanta 
produced during the transition of the geometry from the inflationary to the radiation 
dominated stage of expansion the (canonically normalized) amplitude of fluctuations \& = 
ipa should be defined so that Eq. (|3.17| ) becomes: 



a" 



f + [^ + m V-l]$ = fl. (3.18) 



a 



In the de Sitter stage of expansion Eq. Q3.18|) reduces to 



^' + [A; 2 + = 0, (3.19) 

whereas during the radiation dominated stage of expansion Eq. ( |3.19| ) takes the form 

+ 2 + f j, 2 ( V + 2 Vl ) 2 = Q 

m 
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The solution of Eq. ( [3.19D (with the correct quantum-mechanical normalization for 77 — > 
—00) can be written as 

*i(l) = 4|P H^(-x), (3.21) 



where x = krj and H^> is the first order Hankel function |L3]. In the pure de Sitter case, 



p = 3/2y 1 — (4/9)/i 2 and since p <C 1, p — 3/2; p is a phase factor which has been chosen 
in such a way that 

p = y|e i ^ 1+2 "). (3.22) 



With this choice of p we have that *&i(r)) ~ e ~ lkr ) j y/2k for 77 — ► —00. 

During the radiation dominated stage of expansion Eq. ( |3.20| ) is the equation of 
parabolic cylinder functions |13[]. The solutions turning into positive and negative fre- 
quencies for 77 — > +00 are then 

/r( ^ )= (2^ elfjD -^^ ( * e " 4f2) ' 

= (^T7I e " lf A^|(e- J ^), (3.23) 

where 

z = ^(n + 2m), q = —, (3.24) 

and where £) CT are the parabolic cylinder functions in the Whittaker's notation. The 
solution of Eq. flCTD 

^(77) = c + (% r (r7) + c_(fc)^(?7), (3.25) 

is given in terms of c + (&) and c_(fc) which are the two (complex) Bogoliubov coefficients 
satisfying |c + (/c)| 2 — |c_| 2 = 1. In a second quantized approach |c_(fc)| 2 is the mean 
number of created quanta, whereas in a semi-classical approach c_(/c) can be viewed as 
the coefficient parametrising the mixing between positive and negative frequency modes. 
In the case c_(fc) ~ no mixing takes place and no amplification is produced. In order 
to determine c±(k), ^(77) and *ff T (r)) should be continuously matched in 77 = —rji, namely 

n-m) = (3-26) 

By solving this system, an exact expression for the Bogoliubov coefficients is obtained 
which is, in general a function of two variables : p = 777,771 and X\ = kr\\. Since p <C 1 the 
exact result can be expanded, in this limit, 



c+(k) = ne l *^ ^-^^ S 2 (x 1 ,p)p * + [Si(xi,p) + S 2 (x u p)\p± } + <D(p*), 

c_(k) = ne-^l [ --^^S 2 (x 1 ,p)p-* + ^zVfofap) + S 2 (x 1 ,p))f*} + 0()m) 
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where Si(x\, p) and S^O^i, p) contain the explicit dependence upon the Hankel's functions: 
S 1 (x u p)=e^ 1+ ^H^\x 1 ), 



S 2 (x 1 ,p) = y/x[e l ± 
If p ~ 3/2 Eq. (|3~28j ) gives 

c+(*) 



if(l+2p) 



1 # (1) fa;^ 

2 JX\ 



(3.28) 



c-{k) 



e s V7T 



2T(|)^ 2rg)Aii 



4/^ — V4, 

3 _ 

1 i x l 



i a 
- + 



- l)^/4 1 



2r(fvt 



V^T(j) 
1 (^ + l)/i 1/ 4 i 



si) + ogi 5 m 



In the limit x\ = kr\\ -C 1 the mean number of created quanta can be finally approximated 

as 

n{k) ~ |c_(£;)| 2 = qlkri^p- 1 / 2 (3.30) 

where g is a numerical coefficient of the order of 10~ 2 . The energy density of the created 
(massive) quanta can be estimated from 



dp-4, 



d 3 . 
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m n{k) 



(3.31) 



where u = k/a is the physical momentum. In the case of a de Sitter phase (p = 3/2) the 
typical energy density of the produced fluctuations is 



m 



-1/2 



a 



(3.32) 



The produced massive quanta may become dominant. If they become dominant after 
already decayed they will not lead to further constraints on the scenario. If they become 
dominant prior to the decay of <fi further constraints may be envisaged. The scale at which 
the massive fluctuations become dominant with respect to the radiation background can 
be determined by requiring that p^{rj^) ~ p 7 (^*) implying that 



qmHf 



m 



-1/2 



/Ol 



Hi Ml 



aiV 
a*/ 



which translates into 



2 4 

~ q me , 



(3.33) 



(3.34) 



where e = H^/M-p. In order to make sure that the non-relativistic modes will become 
dominant after already decayed H* < should be imposed, that is to say m > 10 2 
TeV for e ~ 10~ 6 . 
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The maximum tolerable amount of entropy, in order not to wash-out any preexisting 
BAU is model-dependent but, in general, AS < 10 5 seems to be acceptable M, [151, 



Defining as the radiation gas already present at the scale H^, the entropy increase 
from to Td eca y — d 'H^M P is of the order of 

AS=(^^y, (3.35) 

where 

T, = tJ^-) -ro^-i^ (3.36) 



where £ = m/Mp. Demanding that AS < 10 5 implies that 

£ > l(n 10 e 3 . (3.37) 

Taking, as usual, e ~ 10~ 6 , m > 10 GeV. 

Hence, if the constraints pertaining to the homogeneous mode are enforced, the bounds 
coming from the inhomogeneous modes do not invalidate the conclusions of the analysis. 
According to the logic expressed in Eq. ( |3.15| ), an illustrative example is the case where 
m > 10 5 TeV and the BAU is generated after EWPT. In this case the bounds obtained 
in the present section are satisfied and the analysis of the evolution of the inhomoge- 
neous modes shows that the qualitatively new bounds introduced in the picture are less 
constraining than the ones obtained in the analysis of the dynamics of the homogeneous 
mode. 



4 Evolution of the gauge field fluctuations 



The evolution of the field </> during and after the de Sitter stage implies, according to 
Eqs. ( [3.3|) -( |3~9"D , that the two-point function of the gauge field fluctuations may very well 
grow. 

For i] < —rji, a rough approximation suggests that the effect of the ohmic current 
is not present and the gauge field is in the vacuum state with k/2 energy in each of its 
modes. By promoting the classical fields to quantum mechanical operators we have that 
the physical polarizations of the magnetic field can be written as 



d 3 k 



(2vr)3/2 ^ 

where b(krj) = B{krj) / g{rj) obey the equation 



et[a k , a b(kr])e^ s + al k J*(k V )e 



-ik-x] 



(4.1) 



b" + 



k 



f\ 2 



2 - +^ 



0. 



(4.2) 
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Notice that b(kr]) are the correct normal modes whose limit (for 77 — > —00) should be 
normalized to \Jk/2e~ tkri . The two-point correlation function of the magnetic fluctuations 
can then be expressed as 

^(f.ij) = Ux^x + r^)) = J —P^bik^Wik,^, (4.3) 

where 

PnW = ( 4 - 4 ) 

The magnetic energy density, derived from the energy-momentum tensor corresponding 
to the action of Eq. ( |2.2| ), is related to the trace of the correlation function reported in 
Eq. (|4.1|) over the physical polarizations: 

/ \ f / 7 \ sin kv dk .. r . 

pB{r,V) = J PB{k,V)—^ — ^ (4-5) 

where 

P b(M) = ^* 3 |&(M)| 2 - (4-6) 
A necessary condition in order to assess that gauge field fluctuations grow during the de 



Sitter phase is that the two-point function increases in the limit r] — > —771 ||T8| , [19 . 

Suppose that the gauge coupling decreases with monotonic dependence upon the field 
0, namely 



^) = ( 2 1 ^ i ) t , A>0. (4.7) 

This parametrisation is purely phenomenological, however, it allows to take into account, 
at once, some physically interesting cases like the one suggested by the low-energy string 
effective action where, in the limit of 0/Mp < 1, g 2 (4>) ~ 0. 

Using Eq. ( |3.3| ) and Eq. ( |4.7| ) into Eq. ( [4.2| ) the time evolution of the normal modes 
of the magnetic field can be found analytically since the specific form of Eq. (|4.2j ) falls in 
the same category of Eq. ( |3.19| ). Hence, 



b{k,7 ] ) = N^h ] Hl 2 \kr ] ), N = —e- l ^ 1+2u \ (4.8) 

with v = (3A + l)/2 and where H^\kvi) the Hankel function of second kind Notice 
that the normalization iV has been chosen in such a way that for 77 — > — 00 the correct 
quantum mechanical normalization is reproduced. Consequently, following Eq. (|4.3|) , the 
two-point function evolves as 



lim Sij(r,r}) 
r)-+-m 



-3A 

(4.9) 
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Since the two-point function increases magnetic fluctuations are generated. 

For sake of completeness the case of increasing gauge coupling will now be examined 
using the following phenomenological parameterisation 

9(v) = (^ 1 rK S>0. (4.10) 

Again different scenarios can be imagined. For instance, one coould argue in favour of 
scenarios where the gauge coupling depends upon </> (or upon 77) in a highly non-monotonic 
way. For the illustrative purposes of the present investigation it is however sufficient to 
focus the attention on the case of monotonic dependence. 

Following now the same steps outlined in the case of decreasing gauge coupling the 
evolution of the two-point function can be obtained 

2-35 



lim Gij(r,r)) 
n-*-m. 



m 



(4.11) 



lim Gij(r,rj) 



3d 



(4.12) 



for 5 > 1/3 and 

Tj_ 

Vi 

for S < 1/3. If 8 < 1/3 the correlation function decreases and this signals that large scale 
magnetic fields are not produced. 

The back-reaction of the produced fluctuations can be safely neglected in de Sitter 
space. Looking at Eqs. ( f2T3|) and Q2.11|) it can happen that if the magnetic fluctuations 
grow too much the term at the right hand side will become of the same order of the others. 
This is not the case. Using the conventions of this Section together with the explicit form 
of the scale factor in the de Sitter phase it can be shown that 



dg r>2 " - 2 " 



g 3 a 2 



\ki ll \ b - 2u (4.13) 



where v is determined from the specific power dependence of the coupling as a function 
of <p. Since we are interested in large scale modes we have kr\\ <C 1. Therefore the back 
reaction effects are relevant towards the end of the de Sitter phase (i.e. rj ~ — 771 ) and 
for k ~ r/f 1 , namely exactly for the modes not relevant for the present investigation. 

After the end of inflation the onset of the conductivity dominated regime may not be 
instantaneous. In this case after rji the presence of a reheating phase should be taken 
into account. Suppose, for instance, that g{rj) decreases according to Eq. Q4.7|) . Suppose, 
moreover, that during reheating the background is dominated by the coherent oscillations 
of the inflaton. In this case the effective evolution of the geometry will be dominated by 
matter with 0(77) ~ rj 2 . According to Eq. fl3.4p (with a ~ 2), ~ r]~ 3 . The value of the 
magnetic inhomogeneities at r/ r will be given by solving Eq. ((4.2[) 



/> •//, ) ^ -l| ('— "j^ + A 2 f— V~^ A , (4.14) 
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where A\ and A2 are two arbitrary constants. Depending upon the value of A the fastest 
growing solution is selected in Eq. Q4.14 ). This phase may induce further amplification on 
the two-point function since its main effect is to delay the conductivity-dominated regime. 



4.1 Generalized MHD equations 

For rj > i] r the role of the Ohmic diffusion becomes important and the evolution of the 
magnetic inhomogeneities will be described by the MHD equations generalized to the case 
of time varying gauge coupling. 

The ordinary (i.e. fixed coupling) MHD treatment is an effective description valid 
for length scales larger than the Debye radius and for frequencies smaller than the iono- 
acoustic frequency. This means that MHD is accurate in reproducing the spectrum of 
plasma excitations obtained from the full kinetic (Vlasov-Landau) approach but only for 
sufficiently low frequencies and for sufficiently large scales. 

Implicit in the ordinary MHD analysis is the assumption that the plasma has to be 
electrically neutral (V • E — 0) over length scales larger than the Debye radius. Thus, 
this system of equations cannot be applied for distances shorter than the Debye radius 
and for frequencies larger than the plasma frequency where a kinetic description should 
be employed. 

MHD equations can be derived from a microscopic (kinetic) approach and also from 
a macroscopic approach where the displacement current is neglected |Mj . If the displace- 
ment current is neglected the electric field can be expressed using the Ohm law and the 
magnetic diffusivity equation can be derived 

OB 1 

— = V x (v x B) + -V 2 B. (4.15) 
Or] a 

The term containing the bulk velocity field is called dynamo term and it receives con- 
tribution provided parity is globally broken over the physical size of the plasma. In Eq. 
( [4.15] ) the contribution containing the conductivity is usually called magnetic diffusivity 
term. 

In the superconducting (or ideal) approximation the resistivity of the plasma goes 
to zero and the induced (Ohmic) electric field is orthogonal both to the bulk velocity 
of the plasma and to the magnetic field [i. e. E ~ — v X B\. In the real (or resistive) 
approximation the resistivity may be very small but it is always finite and the Ohmic field 
can be expressed as 

E~ X -vxB. (4.16) 

cr 

If the gauge coupling changes with time the system of equations obtained by neglecting 
the displacement current receives new contributions and the relevant equations can be 
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obtained, in the resistive approximation: 



Vx£=--, (4.17) 

E = --$x B, (4.18) 
a 



\vxB = J-2^- 3 E, (4.19) 
g 2 2g s 



Using Eqs. ( |4.17|) -( |4.19| ) the generalized magnetic diffusivity equation can be obtained: 

(1 - 9 -)^— = Vx(vxB) + —,V 2 B. (4.20) 

a g 2 g dr] ag 2 

Notice that Eq. QO0|) reproduces Eq. tfLWj ) if g' -> 0. 

Suppose now that the plasma, whose effective Ohmic description has been presented, 
is relativistic. In the case when the coupling is constant a is constant and it is given by 

<T = (T c (Tj)a(Tl), (4.21) 

where a c ~ T/g 2 scales as the inverse of a(rj) if the evolution of the Universe is, to a good 
approximation, adiabatic. 

If g is not constant, a is not constant anymore but it decreases if the gauge coupling 
increases and, vice versa, it increases if the gauge coupling decreases. In spite of this, in 
the generalized MHD equations, the combination which appears is always ag 2 which is 
roughly constant for an adiabatically expanding Universe. 

In the approximation of instantaneous reheating, the solution of Eq. ( f4.20|) is given 

by 

B(k, rj) = B(k, rji)e ^ 2 ~ 2 ^ . (4.22) 

According to Eqs. ( |3.13| )-( |3.14| ), in order to get the coupling frozen prior to H cw ~ 1CT 17 
GeV, m > 10 5 TeV shall be required. If the gauge coupling is always decreasing as a 
function of r), it can be parametrized by Eq. ( [4.7|) . Hence Eq. (|4.22| ) can be evaluated by 
using the explicit evolution of <fi as obtained from Eqs. (|3.8|)-(|3^) implying that r ~ rf 1 , 
0m ~ ?7~ 3//2 and <p c ~ ri~ 3 . The result is 

B(k, Vo ) =l( Vl , Vm , V ^ Vc )e-^ {m+m) (4.23) 
where t]q is the present time and 

Ita.W,^) = [(^^)(^^)»(^^)r A[ * r (4.24) 
A + ag 2 i] m 3A + 2ag 2 i] c 3A + ag 2 !]^ 
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Concerning Eqs. (|4.23|) -( 4.24 ) few comments are in order. From Eq. (|4.23|) all the modes 



k * >h 2 U JL (4.25) 
Vo 

are suppressed by the effect of the conductivity. The present value of uj a (r)o) [] can be 
estimated by recalling that l/r]o ~ H^a® where Hq ~ 10~ 61 Mp. Thus u) a ~ Hz. 
Present modes of the magnetic fields are dissipated if to > io a . 

As far as the problem of galactic magnetic fields is concerned, the relevant set of scales 
range around the Mpc corresponding to present modes of the magnetic field ujq ~ 10 -14 
Hz, i.e. uj g <ti uj a . 



5 Estimates of large scale magnetic fields 

In the approximation of instantaneous reheating the typical (present) frequency corre- 
sponding to the end of inflation can be computed and it turns out to be 

^M-^T^e 1 / 2 ^- 1 - (5-1) 

Since Td ec ~ 0.26 eV , z^ c T^ cc ~ 100 GHz. Eq. ( phT| ) can be obtained by red-shifting the 
highest mode, i.e. ^i(?7i) ~ Hi through the different stages of the evolution of the model, 
namely, according to Eqs. ( |3.11|) and (|3.14j ), from rji down to f] m and from r] m down to 



r] rmc . Recall that from 7] c to r/^ the Universe is, effectively, matter dominated. The other 
typical frequencies appearing in the time evolution of the gauge coupling can be written, 
in units of u;i(?7o), as 

u m {Vo) _ c -l/2 £—1/2 



e 12 e- 1/2 <p, 

e -l/2 ^7/6 ^2/3^ (5 2) 



where, as in the case of Eq. (|5.1| ) all the frequencies are evaluated at the present time. 

In the case of decreasing gauge coupling [described by Eq. (|4.7|) 1 the amount of 
generated large scale magnetic field can be estimated from Eqs. (|4.2|) - (|4.8|) together with 
Eqs. ( |4.20p — ( 4.23| ). Bearing in mind that the typical frequency scale corresponding to 1 



Mpc is 10 14 Hz we have that the ratio of the magnetic to radiation energy density is: 

tb(oo G ) = /(A) e§ A 10-^-^ T(u G ), (5.3) 



With uj(rj) ~ k/a(r)) the physical momentum will be denoted. 
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where 
and 
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which means, using Eqs. (|5.1|)-(|5T2|). 

T(^ G ) = e^[e- 1 /V- 1 e 5/2 ]" A ^. (5.6) 

Notice that in Eq. ( |5.5| ) is the present CMB temperature. 

In order to illustrate the regions of the parameter space where magnetogenesis is 
possible cp ~ 1 will be assumed. In Fig. [I] and Fig. |2| the case of decreasing gauge 
coupling is discussed. In Fig. p], A is fixed and the exclusion plot is given in terms of e 
and £. In particular, for illustration, A = 1 has been chosen. The choice of A ~ 1 implies 
that g 2 (4>) ~ <f). Such a case is favoured from the tree- level string effective action where 
the effective coupling can be approximated as g 2 (4>) ~ 0/Mp for < Mp. 

The two vertical lines lines mark, respectively, the bounds coming from BBN [i.e. 
£ > 10~ 15 , from Eq. (|3.14|) 1 and from the electroweak epoch [i.e. £ > 10~ n , from Eq. 



( p.!5| )]. For consistency with the assumptions of Eqs. (|3.2|)-(JP|) m/ Hi <C 1, implying 
e ^> £. With the lower dashed line the bound e > £ is reported. Notice that the 
requirement of Eq. ( |3.37| ) is not numerically relevant. The upper dashed line is obtained 
by requiring, according to Eq. ( |3.34| ), that > H*. Finally, the full (diagonal) line is 
derived by imposing on Eqs. (|5.3|) - (|5.6|) the dynamo requirement. 

In Fig. the value of £ has been fixed to 10 -11 , as required by the considerations 
related to the electroweak epoch and the magnetogenesis region is described in terms of 
e and A, both varying over their physical range. The two horizontal lines fix the bounds 
coming from e < and from e > £. With the dashed line the curve tb{ojg) ~ 10 -20 is 
denoted. Hence, values larger than the dynamo requirement are allowed. This observation 
may be relevant in the context of magnetic fields associated with clusters. In Fig. ^ A lies 
in the range < A < 4/3. This choice guarantees the growth of the correlation function 
of the magnetic inhomogeneities during the de Sitter stage. 

In order not to conflict with large scale bounds coming from the isotropy of the CMB 
the energy spectra of the produced gauge field fluctuations have to decay at large distance 
scales, implying that < A < 4/3. To be compatible with CMB anisotropies r B (o;dec) < 
1CT 10 should be imposed (where oj^ec — 1CT 16 Hz). If the condition < A < 4/3 is 
enforced, the spectra increase with frequency and the possible bounds coming from the 
anisotropy of the CMB are satisfied. 

An analogous estimate can be obtained in the case of increasing gauge coupling dis- 
cussed in Eq. ( |4.1U| ). In this case 

tbM ~ mel 6 - 1 V 4 - 25 10- 25{6 - 35) T(uj g ), (5.7) 
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Example for fixed X 




Figure 1: The shaded area illustrates the region where magnetogenesis is possible in the 
case where A = 1 and tp ~ 1. The vertical lines correspond to the requirements coming 
from BBN and from the EWPT epoch. 
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Example for fixed § 
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Example for fixed 8 
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Figure 3: The magnetogenesis region in the case 5 = 2 (i.e. increasing gauge coupling) 
and tp ~ 1. The exclusion plot is then given in terms of e and £. 



where 



and 



/(*) 



o3<5-3 / o i 



7T" 



.2 2, 

r(^ G ) = e ^[ e - 1 /V 1 e 5/2 ]^. 



(5.8) 
(5.9) 



In Fig. |^ and Fig. |] the requirements coming from the dynamo mechanism as well as the 
other theoretical constraints are illustrated. For both plots e < 10~ 6 and <p ~ 1. In Fig. 
H the case 6 = 2 is illustrated. The shaded area selects the region of the parameter space 
where the dynamo requirement is imposed on Eq. (|5.7| ). The two vertical lines illustrate 
the conditions of Eqs. (|3.14|) -( |3.15|) . As in Fig. |2| the two dashed lines correspond to the 
constraints coming from the inhomogeneous modes and from the condition e > £. 

The requirement that the spectra decrease at large distance scales implies, in this case 
that 5 < 2. The condition on the growth of the correlation function obtained in Eqs. 
( |4.11| )- (|4.12| ) imply 5 > 1/3. Thus the interesting physical range of Eq. fl4.10|) and (|5.7|) 
will be 1/3 < 5 < 2. 

In Fig. II the parameter space is illustrated for fixed values of £, i.e. £ ~ 1CT 11 . As in 
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Example for fixed ^ 




Figure 4: The magnetogenesis region is illustrated in the case of increasing gauge coupling 
in the (e, 8) plane. Notice that in the present example £ ~ 1CT 11 (for compatibility with 
the EW epoch) and ip ~ 1. 

Fig. |3] the full and dashed curves correspond to the dynamo requirements imposed on Eq. 
(p.7[). The shaded area selects the allowed region in the space of the parameters where 
magnetogenesis is possible for 10~ n < e < 10~ 6 . As in the case of Fig. |2| there are regions 
in the shaded area where values much larger than the magnetogenesis requirement are 
possible (see the dashed line in Fig. |]). 

As it has been pointed out in deriving the theoretical bounds on the scenario the 
requirement £ > 10~ n may be too restrictive since it excludes the variation of the gauge 
coupling at the electroweak time. To relax this assumption is possible and it would 
require a precise analysis of the dynamics of the EWPT in the presence of time varying 
gauge coupling. At the moment this kind of analysis is not available. Summarizing this 
illustrative discussion, there are regions in the parameter space of the models wher all the 
theoretical constraints are satisfied and where magnetogenesis is possible. In particular, 
it is possible that the gauge coupling freezes prior to the electroweak epoch, leading still 
to magnetogenesis. 
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6 Concluding remarks 



There are no reasons why the gauge couplings should be constant throughout all the 
history of the Universe. If they are allowed to change prior to the formation of the light 
elements they can lead to computable differences in the cosmological evolution. 

In the present paper the interplay between inflationary magnetogenesis and the evo- 
lution of the Abelian gauge coupling has been addressed. In a phenomenologically rea- 
sonable model of inflationary and post-inflationary evolution the relaxation of the gauge 
coupling leads to a growth in the correlation function of magnetic inhomogeneities. Large 
scale magnetic fields are then generated. The evolution of the gauge coupling is driven, 
in the present context, by a massive scalar. 

Since the gauge coupling evolves significant changes in the evolution of the plasma can 
be envisaged. In the present investigation the ordinary MHD equations have been gener- 
alized to the case of time evolving gauge coupling but other effects could be envisaged. 
In particular, the generlization of the full kinetic approach to the case of time evolving 
"electron" charge would be of related interest. 

The value of large scale magnetic fields produced with this mechanism has been es- 
timated. For a broad range of parameters the obtained values of the magnetic fields are 
much larger than the dynamo requirements. 

In the present investigation the main assumption has been that the only gauge cou- 
pling free to evolve is the Abelian one. Furthermore, the parameters of the model have 
been chosen in such a way that the gauge coupling is not dynamical by the onset of the 
electroweak phase transition. It would be interesting to relax both assumptions since they 
may lead to potential differences with the standard scenarios. 

Therefore, in many respects, the present investigation is not conclusive. At the same 
time it shows that acceptable models for the evolution of the gauge couplings can be 
obtained in a standard cosmological framework. 
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